Dimerization versus Orbital Moment Ordering in the Mott insulator YVO3 
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We use exact diagonalization combined with mean-field theory to investigate the phase diagram of 
the spin-orbital model for cubic vanadates. The spin-orbit coupling competes with Hund's exchange 
and triggers a novel phase, with the ordering of t2 g orbital magnetic moments stabilized by the tilting 
of V06 octahedra. It explains qualitatively spin canting and reduction of magnetization observed 
in YVO3. At finite temperature an orbital Peierls instability in the C-type antiferromagnetic phase 
induces modulation of magnetic exchange constants even in absence of lattice distortions. The 
calculated spin structure factor shows a magnon splitting due to the orbital Peierls dimerization. 
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Many transition metal oxides are Mott-Hubbard in- 
sulators, in which local Coulomb interaction oc U sup- 
presses charge fluctuations and leads to strongly corre- 
lated 3d electrons at transition metal ions |l| . When de- 
generate d orbitals are partly filled, the orbital degrees of 
freedom have to be considered on equal footing with elec- 
tron spins and the magnetic properties of undoped com- 
pounds are described by spin-orbital superexchange (SE) 
models 0, 0] • Such SE interactions are typically strongly 
frustrated on a perovskite lattice, leading to enhanced 
quantum effects ;4]. In systems with e g orbital degener- 
acy (manganites, cuprates) this frustration is usually re- 
moved by a structural transition that occurs well above 
the magnetic ordering temperature and lifts the orbital 
degeneracy via the cooperative Jahn- Teller (JT) effect. 

A different situation arises when ti g orbitals are partly 
filled like in titanium and vanadium oxides. As the JT 
coupling is much weaker, the intrinsic frustration be- 
tween spin and orbital degrees of freedom may show up in 
this case. Unusual magnetic properties of titanates have 
recently been discussed in terms of coupled spin-orbital 
SE dynamics 0| • In addition to SE, spin and orbital occu- 
pancies of tig levels are coupled also via atomic spin-orbit 
interaction, H so oc A(S ■ I), which is particularly relevant 
for vanadates with a triplet 3 T 2 ground state of V 3+ ions. 
As a result of intersite SE and on-site A interactions spin 
and orbital orderings/fluctuations are strongly coupled, 
as observed in the canonical spin- orbital system V2O3 
Q, as well as in cubic LaVC-3 Q. 

The magnetic properties of YVO3 are particularly puz- 
zling |8|], and indicate dimerization along the FM direc- 
tion within the C-AF phase [jj]. In this Letter we argue 
that such an exotic C-AF phase follows from the realis- 
tic spin-orbital model for vanadates that emphasizes the 
competition between SE bond physics and intraatomic 
spin-orbit coupling cx A. We investigate the phase dia- 
gram of this model and show that orbital moments are 
induced in the C-phase by finite A, and form at larger A 
a novel orbital moment (OM) phase. 



The superexchange in cubic vanadates originates from 
virtual charge excitations, dfd 2 — ► rffdj, by the hopping t 
which couples pairs of identical orbitals. When such pro- 
cesses are analyzed on individual bonds (ij) | 7 along 
each cubic axis 7 = a, b, c, one finds the spin-orbital 
Hamiltonian with S = 1 spins (J = At 2 /U) [l(|, 
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where the orbital operators J^P and depend on the 
pseudospin r = 1/2 operators = {rf ,rf ,rf}, given 
by two orbital flavors active along a given direction 7. 
For instance, yz and zx orbitals are active along c axis, 
and we label them as a and b, as they lie in the planes 
orthogonal to these axes. The general form of the su- 
perexchange Q was discussed before, and we have shown 
that strong quantum fluctuations in the orbital sector 
provide a mechanism for the C-AF phase ^(J. When c 
(xy) orbitals are occupied (n, c = 1), as suggested by the 
electronic structure [ll( and by the lattice distortions in 
YVO3 t the electron densities in a and b orbitals 

satisfy the local constraint ni a -\-nn = 1. The interactions 
along the c axis simplify then to: 
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they involve the fluctuations of a and b orbitals oc ?i ■ tj , 
and fi<SiTj — T?Tj — rfrj + t'tj, while the interactions 
along the 7 = a(b) axis depend on the static correlations 
cc nibUjb (rii a nja) only; for instance: 
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The Hund's exchange 77 = Jh/U determines the mul- 
tiple! structure of d 3 excited states which enters via 
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the coefficients: R = r]/(l — 3r]) and r = 77/(1 + 27/). 
The pseudospin operators in Eqs. ©-I© may be rep- 
resented by Schwinger bosons: rf = {a\b i + feJaJ/2, 
t? = i{a\bi - ftJcO/2, 7f = (n m - n J/2. 

The individual V06 octahedra are tilted by angle fa = 
±0, which alternate along the c axis As the icy 

orbital is inactive, two components of the orbital moment 
li are quenched, while the third one (2rf), parallel to 
the local axis of a VOg octahedron, couples to the spin 
projection. Because of AF correlations of rf moments, 
spin-orbit coupling induces a staggered spin component. 
As the spin interactions are FM, weak spin-orbit coupling 
would give no energy gain, if the spins were oriented along 
the c axis. Thus, finite A breaks the SU(2) symmetry and 
favors easy magnetization axis within the (a, b) plane. As 
quantization axis for li (Si) we use the octahedral axis 
(and its projection on the (a, b) plane), respectively. The 
spin-orbit coupling in Eq. JJJ is then given by: 

H so = 2A^2 (Sf cos fa + Sf sin fa) rf , (6) 

i 

and we use A = A/J as a free parameter. In order to 
understand the important consequences of the tilting for 
the interplay between spin and orbital degrees of freedom, 
we consider first the idealized structure with = 0. The 
coherent spin-and-orbital fluctuations lower then the en- 
ergy due to on-site correlations (Sfrf) < 0. Since these 
fluctuations do not couple to the spin order, no orbital 
moments can be induced at small A as long as <j> = 0. 

Even in the absence of the spin-orbit term (A = 0), 
the vanadate spin-orbital model poses a highly non- 
trivial quantum problem. We obtained first qualitative 
insight into the possible types of magnetic and orbital 
ordering by investigating the stability of different phases 
within the mean-field approximation (MFA), but includ- 
ing the leading orbital fluctuations on FM bonds along 
the c axis. In the absence of Hund's exchange (ry = 0), 
two orbital flavors experience an antiferro-orbital (AO) 
coupling on these bonds due to jff Eq. J5J, but are de- 
coupled within the (a, b) planes (K^^ — 13] ) . This 
one-dimensional (ID) system is unstable towards dimer- 
ization with orbital singlets and FM interactions at every 
second bond along the c axis ^J, stabilizing the orbital 
valence bond (OVB) phase. The spin interactions J^'^ 
Eq. Q and the intersinglet interactions along c axis are 
AF. In contrast, for large 77 more energy is gained when 
the orbital singlets resonate along the c direction, giving 
uniform FM interactions in the C- AF phase [l0| . We de- 
termined the quantum energy due to orbital fluctuations 
using the orbital waves found in the Schulz approxima- 
tion, known to be accurate for weakly coupled AF chains 
[lfij . Using this approach, the transition from the OVB 
to C-AF phase [(Sf) = S z e lR ^ c with Q c = (tt.tTjO)] 
takes place at r\ c ~ 0.09 (Fig. [TJ, and the orbital or- 




FIG. 1: Phase diagram of the spin-orbital model in the (77, A) 
plane at T = 0, reflecting the competition between orbital 
valence bond (OVB), staggered orbitals (C-AF) and orbital 
moment ordering (OM), as obtained by the ED of a four- 
site embedded chain for <f> — 11° (circles), and in the MFA 
(dashed lines). Orbital moments in the OM phase (violet ar- 
rows) induce spin canting (blue arrows) with angle ip. (YVO3: 
r; ~ 0.12, A ~ 0.3 - 0.4 [ill). 

dering, (rf) = T z e lR iQG w ith Qc — (tt,it, it), sets in, 
promoted by the AO interactions Eq. (J5J. 

The ground state changes qualitatively at finite A and 
4> > 0. The magnetic moments (Sf) induce then the 
orbital moments, 2(rf) = l z e lR iQA^ which stagger along 
the c axis with Qa — (0, 0, tt). This novel type of ordering 
with If ^ can be described as a staggered ordering 
of complex orbitals a ± ib (corresponding to If = =pl 
eigenstates) that competes at 77 > r\ c with the staggered 
(a/b) ordering of real orbitals with t z ^ in the C-AF 
phase. Already at small A the orbital moments induce in 
turn opposite to them weak (Sf) ^ moments, lowering 
the energy by (Sfrf) < 0. In the OM phase favored 
at large A (Fig. [TJ, the spin order has therefore two 
components: (Sf) = S z e iS *® a , and (Sf) = S x e iRi ^ G . 

An unbiased information about the spin and orbital de- 
grees of freedom was obtained by the accurate treatment 
of quantum effects within the exact diagonalization (ED) 
method. Thereby the coupled spin-orbital excitations, 
terms oc SfS^T^T^ in Eq. ijljl. are now fully included. 
We performed ED of four-site chains along c axis, both 
for free and periodic boundary conditions (PBC). 

We were surprised to see that the exact ground state 
of a free chain at r\ = A = consists indeed in a very 
good approximation of two orbital singlets on the ex- 
ternal (12) and (34) FM bonds {(n ■ f i+ i) = -0.729), 
connected by an AF bond (23) with decoupled orbitals 
({n ■ T i+ i) = -0.038). The spin correlations are FM/AF 
on the external/central bond, (Si ■ Sj+i) = 0.95 and 
— 1.56. With increasing 77 the AF interaction weakens, 
the sequence of spin multiplets labelled by the total spin 
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FIG. 2: Pair correlations along c axis: (a) orbital (t; • 
(b) spin (5i • Si+i), and (c) spin components (S z ) (full lines) 
and (Sf } (dashed, long-dashed lines), as functions of A for the 
OVB (A < A c ~ 0.37) and the OM (A > A c ) phase, found by 
the ED at 77 = 0.07. FM/AF bonds (ij) in the OVB phase 
are shown by solid/dashed lines in (a) and (b). 



St is inverted at r\ c ~ 0.12, and the ground state changes 
from a singlet (St = 0) to a high-spin (S t — 4) state. 
At finite A no level crossing occurs, but the nondegener- 
ate ground state describes a smooth crossover in the spin 
and orbital correlations with increasing 77. We verified 
that several excited states lie within 0.1J away from the 
ground state — all of them would contribute to thermal 
fluctuations already at temperatures T ~ 30 K. 

We simulate a cubic system by including infinitesimal 
symmetry-breaking dimerization field which favors the 
orbital singlets on the bonds (12) and (34) in a cluster 
with the PBC, embedded within one of three phases sta- 
ble in the MFA (Fig. [IJ, with mean-fields determined 
self-consistently in each phase. All phases are character- 
ized by finite magnetic moments (Sf), either staggered 
pairwise (OVB phase), or aligned (C-AF phase) along c 
axis, and weak (Sf) moments. In addition, the orbital 
ordering (l z ^ 0) appears in the OM phase, while the 
orbitals stagger (t z ^ 0) in the C-AF phase. 

By computing the energies of different phases we ob- 
tained the phase diagram that confirms the qualitative 
picture extracted from the MFA (Fig. QJ. All transitions 
are accompanied by reorientation of spins (Figs. [21 and 
[3J. The orbital and spin fluctuations change only weakly 
at small values of A, weak (Sf) moments are ordered 
pairwise on the bonds (12) and (34), and the spin corre- 
lations on the intersinglet bonds are almost classical (Fig. 
[2J . These features show that the spins and orbitals are al- 
most decoupled and the OVB phase is robust. At A > A c 
the on-site correlations (Sfrf) < dominate, while the 
orbital fluctuations are suppressed, and the correlation 
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FIG. 3: Crossover from the C-AF to the OM phase for in- 
creasing A at 77 = 0.12: (a) orbital order parameters: r z 
(dashed line), and V (solid lines) at cj) = 11°; (b) spin canting 
angle ij) — arctan(5 ,I /5 z ) (see Fig. 0, as obtained for the 
tilting: (j> — 5°, 11° and 20° (dashed, solid and dotted line). 



functions (fj • f^+i) approach the classical value —i. As 
the staggered spin moments (Sf) are induced, the spin 
correlations (Si ■ Si+i) become soon AF within the OM 
phase. In this regime the spins follow the spin-orbit cou- 
pling A, and the FM interaction J c is frustrated. 

At the transition from the C-AF to the OM phase 
the orbital ordering changes from staggered real orbitals 
(t z ^ 0) to staggered orbital moments (l z ^ 0), as shown 
in Fig. EJa). As a precursor effect of the forthcoming 
OM phase, the orbital moments If are induced already 
in the C-AF phase by increasing A. The transition to 
the OM phase is accelerated by the increasing tilting an- 
gle <j> [Fig. Ofb)]. Also the spin correlations change here 
discontinuously at the transition (not shown), similar to 
the OVB/OM transition [Fig. Gib)]. The staggered spin 
components in the OM phase (Sf) are similarly large to 
those shown in Fig. Efc) for smaller 77, and the spin cant- 
ing angle if) approaches -| — 4> in the regime A> 1. For re- 
alistic parameters for YV0 3 : J - 30 meV ;10], 77 - 0.12 
(estimated with Jh = 0.64 eV and intraorbital element 
U = 5.5 eV for V 2 + ions 0), and A ~ 0.3 - 0.4 (con- 
sidering A ~ 13 meV for free V 3+ ions 0|)> one finds a 
competition between the staggered a/b orbital order pa- 
rameter (r z ~ 0.25) and the orbital magnetic moments 
(l z ~ 0.30 — 0.35). This reflects the interplay between 
intersite SE and on-site spin-orbit couplings, and hence 
orbital and spin moments are not collinear (except for 
large A values), — in contrast to the conventional pic- 



ture where orbital moments induced by A coupling arc 
antiparallel to spin, as suggested e.g. for V2O3 [19j . 

Finally, we turn to finite temperatures. While the C- 
phase cannot dimerize at T = 0, it dimcrizes at finite T 
due to the intrinsic instability towards alternating orbital 
singlets The spin correlations (Si ■ Si+i) [Fig. E^a)], 
found using open boundary conditions, alternate between 
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FIG. 4: Dimerization in the C-AF phase at T > 0: (a) spin- 
spin correlations (Si ■ Si+i) on strong and weak FM bonds 
(solid and dashed line); (b) spin response S(q, uj) in the dimer- 
ized C-AF phase for q = 0,7r/2,7r. Inset: filled (open) circles 
indicate strong (weak) features in S(<\, uj); lines show the fit- 
ted spin- wave dispersion. Parameters: r\ = 0.12, X = 0.4. 

strong and weak FM bonds due to the orbital Peierls 
dimerization, 2S T = \{fi ■ t^+i) — (fj ■ i ) | , which has a 
distinct maximum at T ~ 0.24 J for r\ = 0.12. Consistent 
with our discussion above, the modulation of (Si ■ Sj+i) 
vanishes in the C-phase at low T. 

Up to now, the sole experimental evidence for dimer- 
ization of the C-phase is the splitting of FM spin waves 
in the neutron scattering study of Ulrich et al. 0. Fig. 
4(b ) shows the dynamical spin structure factor S(q,u>) 
[20j obtained by exact diagonalization of a 4-site cluster 
with PBC at T = 0, assuming the same orbital dimer- 
ization (ri ■ Ti+i) as found above for T/J — 0.25. At 
q = 7r/2 we observe a splitting of the spin wave simi- 
lar to experiment. The finite energy of the q = mode 
results from the A-term and the mean-field coupling to 
neighbor chains. Additional features seen in S(q,u>), e.g. 
for q = it at uj ~ 1.25J, we attribute to the coupling to 
orbital excitations. The spin-wave energies can be fitted 
by a simple Heisenberg model with two FM coupling con- 
stants: J c i — 5.7, J C 2 = 3.3 meV, and a small anisotropy 
term, as shown in the inset (solid lines). Although these 
values are strongly reduced by spin-orbit coupling Eq. 
|JB|, they are still larger than those extracted from the 
spin waves in YVO3: J™ p = 4.0 and J^ p = 2.2 meV at 
T = 85 K We attribute this overestimate of exchange 
interactions to quantum fluctuations involving the occu- 
pancy of xy-orbitals; this will be treated elsewhere. 



Summarizing, we have shown that the spin-orbit cou- 
pling A competes with Hund's exchange in the spin- 
orbital model for cubic vanadates. It leads to a new or- 
bital moment ordered phase at large A and can explain 
qualitatively the spin canting and large reduction of mag- 
netization in the C-phase at smaller A. We argue that 
the ID orbital Peierls instability observed recently in the 
C-AF phase of YVO3 [2( (along the c axis) emerges from 
a combination of quantum effects due to orbital moments 
at A > 0, and thermal fluctuations which favor dimerized 
orbital and spin correlations. 
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